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APPRO'"' '^  POR  P'JBLIC  RJUlSt 

aisiRs.B’jno'j  uirLijuia) 


laJLTlPATH  IN  THE  THREE-DIMENSIONAL  UNDERWATER  ARiLW 

A  planar  array  of  widely  dispersed  hydroprioiicy  doploye>d  in  a  horizontal 
plane  several  hundred  meters  below  th.''  surface  of  the  sea  -been  analyzed, 

and  reported  earlier  [1,  2].  Becausv^  of  the  dispersive  nature  of  the  under¬ 
water  medium,  signal  energy  approach.es  the  array  from  a  source  along  a  n  umber  of 
refracting  paths.  This  results  in  n-.u.ltipath  inter ference  a'^ci  a  possible 
consequent  loss  of  output.  Ray  arrivals  from  th.e  source  typically  fall  into 
a  range  of  vertical  angle.s  at  the  receiver  which  are  t  10°  relative-*  to  tlie 
liorizontal.  The  vertical  beamwidth  of  a  planar  a.ilcnna  is  large*  enough  to 
accept  all  rays  in  sucli  a  range,  heuC'e  it  is  multipath  sensitive.  A  three- 
dimensional  array  is  capad>le  of  a  sh:*>'p.,?r  vt*rLicci].  foi  uf,  and  will  be  lo.us  sensi¬ 
tive  to  the*  kind  of  multipath  typical  in  Ltiis  appl  i  ca  t  iot> .  in  fart,  by  suitably 
processing  the  array  out[iut  there  is  the  p  »sr- i!)i  1  i  ty  that  tlie  mtiltipath  arrivals 
can  be  separate}'/  received  and  then  eonbim^d  in  yhnr.o  t(»  achieve  an  ’’angle  cC 
arrival  diversity”  system^as  suggrjstt^d  in  Figure  3.1.  Such  syster^s  have  been 
proposed  for  tropospheric  su  alter  rec,eivcrs.  - 

We  analyze  the  mean  properties  ''f  suc!i  a  svi-tv-m  below.  If  tlu*  array  is 
focused  to  look  i:\  the  y-z  plane  its  resiyon  n*  to  a  signal  at  angle  y  relative 
to  the  x-z  plain.'  is 
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FIGURE  3.1  IIULTIPLE  VERTICAL  BE*AMS  }'OR  MULTIPATH  RESOLUTION 
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Here  0  is  the  vertical  angle  to  which  the  array  is  to  be  focused; 


(x  ,  y  ,  z  )  is  the  position  of  the  nth  array  element,  (1)  assumes  M  ray 
n  n  n  ' 


arrivals  each  given  by  B  e  with  vertical  arrival  angle  0  . 

m  m 


We  concentrate  on  tije  output  when  llie  source  is  on  the  main  beam;  tliat 


is,  when  (J>  =  90®.  Then 


A(?.0  )  =  I  ):  B, 
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n  m 


m-1  n=l 


QPR  No.  31 


We  calculate  the  mean  power  response  of  the  array  given  by 
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the  Kronecker  delta,  so  that 


=  ",  f,  ?,  <v><«  ■" 

m-l  n^=i  n2=l  ^ 


jk(y  -y  )(sinO  -sinO  ) 
n  T  n  2  ^  s 


jk(z  ”7  )  (cos9  -cosO  ) 

n^  n2  ra  s 


> 


M  .  . 


M  M  jk  y  (sinO  -  sinG  ) 
N  N  .  ^n,  '  m  s 


m=l  [_  n^  =  l  n,  =  l  ^ 


^3^  ^2  ji'-  7.  (cosO  -  cosO  ) 

in  ®  \ 

*  / 


“iV:v  (sinf-  -  sinO  ) 


-ikz  (cos^-  -  cosO  ) 
1^2  m  s 


(5) 


Ol’U  ■)] 


We  assume  all  vectors  (v  ,  z  ) ,  n .  =  1,  2 ,  . . .N  identically  distributed, 

'  n .  n .  1 

1  1 

Furthermore,  the  random  variables  y  and  z  are  assumed  independent  and 

n .  n , 

1  i 

symmetrical  around  the  origin.  Then 


m=l 


N  +  (N  -  N)  • 


<iky  (sin0  --  sin9  )  >. 
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The  expectations  on  the  right  are  characteristic  functions. 


«(jt)  =  ^  (7) 

where  the  random  variable  u  is  cither  v  or  z  and  t  is  correspondingly 

^  n  n  r  o 

either  k(sinO  -  sinO  )  or  k(cosO  -  cosO  ).  The*  random  variables  will  here 
m  s  ^  ra  s 

be  specified  as,  either,  uniformly  distributed  in  an  interval  (-h,  h)  ,  or 

2 

normally  distributed  around  zero  with  variance  O  .  Thus  for  the  uniform 
case 


.  ,  .  sin  ht 

Mjt)  = 


and  for  the  normal  case 
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Jf  the  variables  y  and  z  are  both  ntunial  with  variance  D,,  and 

n  n  y  2 

respect  i  v(' 1  y ,  we  havi' 
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-[a  ^k^(sin9  -sinO  )^  +  o  ^k^(cos9  -cos0  )^] 
y  m  s  ^  ^  q'J 
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If  the  variable  y  is  normal  with  variance  a  and  the  variable  z  is 
n  y  n 

uniform  in  (-*h,  h)  then 


-cr  ^k^(sinO  -sinO 

y  ms 
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An  inspection  of  (10)  or  (11)  leads  to  the  conclusion  that  if  the 
vertical  dimension  of  the  airay  is  in  the  order  of  10  wavelengths  the 
vertical  be^a'.'widtl'i  will  be  about  1°.  r\irt  le.rniore  rays  entering  through 
this  narrow  bea.r;wicth  (hence  oxclud  ie.g  c  lh^r  rays  arr  iving  at  vertical 
angles  outside  the  vert  i  cal  benmwidtli)  wil  1  t>e  suf  f  ic  lent  ly  compact  to 
avoid  the  effect  of  phar>e  decorrel  it  ion  across  the  array. 

Wc  arc  therefore  led  to  propose  the  following  concept.  Let  the 
array  simul tancc'usly  form  contip.uous  vertical  bean5;  as  indicated  in 
Figure  3.1.  Outputs  cor  respond  ing  to  e:ich.  beam  will  be  simultaneously 
present.  Ihese  outputs  are  then  coluTently  combined.  The  operations 
required  are  as  indicated  in  i'igure  3.2.  The  mechanism  being  suggested 
is  similar  to  that  used  in  ang.lo  of  arrival  diversity  coirmuinication  sys¬ 
tems  w'ith  maximal  rat  it  j  c  ordain  ing  of  the  diversity  ignals. .  As  a  rule 
in  these  svstt'ins  each  diversity  l-ranch  has  a  separ.»te  directive  sensor 
and  [)  rtari'^  1  i  i’ i  e  r .  Here  sens(.»rs  and  pri'Ctr:plit  it'rs  are  con'mon  for  all 
branefK‘s. 
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Because  sensors  are  common  one  may  question  the  effect  of  noise 
generated  at  the  sensor  or  preamplifier  input.  Will  such  noise  be 
independent  when  observed  at  the  point  of  combination  of  the  diversity 
branches?  The  following  is  a  discussion  of  that  puint- 
A  filter 


H(f)  -  A(f) 

which  acts  as  a  constant  gain  device  and  constant  phase  shifter  -  that  is, 
with 


A(f)  =  A 

(()(f)  =  “(t  sgn  f  a  constant) 

can  be  represented  by 

H(f)  =  A  e  ^  =  A  cos  (-:  sgn  f)  +  j  A  sin(-v  sgn  f) 

=  A  cos:  - j  A  s in  :  s  gn  f 

A  filter  with  f r eq no n c  y  c  hn  r a c  t  o r i s  t i c 

fl^^(f)  =  -j  ^gn  f 

is  a  Hilbert  transforming  filter  so  that  a  wave  function  n(t)  applied  to 
H(f)  as  defined  abo\^e  emerges  as 

n  (t)  =  A  cosh  n(t)  -f  A  sinh  n(t) 
o 

where  n ( t )  is  the  Hilbert  Transform  of  n ( t ) . 

By  dirret  .Triplication  of  the  definitions  and  by  use  of  the  statistical 
properties;  of  the  Hilbert  irar.sfnr’u  one  can  s-how  that  lor  stationarv^ 
zero  niean ,  j>rocessc's, 


QPK  No.  31 


^n^(t+T)  n^(t)^  =  A^<^n(t+T)  n(t)^  = 

Input  and  output  autocorrelation  functions  are  proportional. 
the  input  autocorrelation  function.  Also, 

^  n  (t+T)  n(T)  \  A  cos(P  R  (t)  +  A  s±n(p  R  (t)  , 
x  o  /  n  n 

giving  a  relationship  between  input-output  cross-correlation  function  and 

the  input  autocorrelation  function.  R  (t)  is  the  Hilbert  Transform  of 

n 

Rn(t).  Finally  for  an  input  n(t)  applied  to  separate  filters 


Hi(f)  =  e 


sgn  f 


and  H^Cf)  =  A^  e 


-3^2  ^ 


the  cro.ss-corrclat  i  3n  function  of  the  two  outputs  n  (t)  and  n  (t)  is 

Oi 


/  n  (t+T)  n  (t)N  =  A  A  [cos(o 


-4)^)  R^(T)  -f 


Consider  now  the  block  diagram  of  Figure  3.2.  Each  branch  is  com¬ 
prised  of  the  sum  of  ll  inputs,  one  from  each  array  element  and  phase 
shifter.  Brancii  1,  for  instance,  contains  a  wave  function 


N  (t)  =  I  n  ,(t)  =  I  [cosi  11  (t)  +  sinA  n  (t)] 
^  n=l  n=l  " 


whore  n  (i)  is  the  output  at  the  n’th  sensor,  and  n  ,  (t)  is  the  phase- 
n  onl 


shifted  output  of  the  n’th  sensor  contributing  !:o  branch  1.  is  defined 

in  Kip^urc  3.2. 
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The  output  N^(t)  will  ultimately  be  applied  through  the  second  filter 
with  characteristic  given  by  resulting  in  an  output 


N 

M^(t)  ^  I  [cos(^^^^  +  sin((Jj^^+  li^^) 

n=l 


For  all  I  branches  taken  together  we  get 

I  IN 

N(t)  =  I  M  (t)  ^  I  I  A  [cos(ij) 

^r=i  ^  i=l  n=l  ^  ^ 

+  sin(;^j.+  n^(c)] 

/  i 

We  view  the  n  (l)  to  be*  Gaussian  noise,  present  at  the  sensor  outputs 
n 

generated  in  the  sensor,  its  I'r  ear.;)]  i  f  j  or  an,d  trie  sensor  *s  iimaedicate 

surroundings.  V.’e  a.'^suire  n  (t)  is  indepcad'nt  of  n  (t)  for  n  ^  in.  Thus 

n  I'll 

we  exclude  external  noise  which  my  bv-  corri'lated  across  several  sensors 
The  variance  of  h'(t)  is  then 
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Wc  have  in.ed  the  i  lol  -  youd-  ;u  •  Mdtf  i>  ..  .r,  :  ,})  .r  i  :id»  •juaidencc* 


*,o  . 


of  n  (t)  and  n  (t)  for  all  n  and  m  for  a  Gaussian  process.  It  can  be 
n  tn 

shown  that 


SO  that 


<(N^t))>  =  <(a.A.  cos(^^.-  4>^)y 


To  continue  this  analysis  ve  require  the  joint  statistical  properties 
of  the  phase  shifts  and  the  anplitude  factors.  For  our  purposes  at  present 
we  may  assume  the  constant  for  all  i  1,  2,  ...I.  But  the  phase  shift 
properties  are  needed.  h’ote  the  following.  If  =  1,  all  i,  and 

<;cos(9ni-  =  0. 

except  when  i  =  j  >  then 

<N^t)>  -  I  5 

However,  if  the  angular  differences  were  small  so  that 


=  1 


for  al 1  j  and  j  then 

<b'(0  >  ■  I  <%'(‘)> 

In  t iiv^  latter  thit  branch  noises  »are  correlated  and  add  coherent Iv, 

In  tlu*  forner  case  tlicy  are  unco v related  and  add  incoherently . 

The  difference  angle  statistics  are  under  investigation  and  will  be 
reported  later.  A  preliminary  calculation  indicates  that  for  the  array 
size  onv  I  :>  j  i)ncd  in  thi.s  application  the  invyilar  diffc’rcncos  may  be  large 
enough  rc>!  t!ie  tirst.  crndirih'jn  above  to  be  aj'tprox  imal  el  y  correct. 
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